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We present a theoretical study of nonlinear pattern formation in parametric surface 
waves for fluids of low viscosity, and in the limit of large aspect ratio. The analysis is 
based on a quasi-potential approximation to the equations governing fluid motion, fol- 
lowed by a multiscale asymptotic expansion in the distance away from threshold. Close 
to onset, the asymptotic expansion yields an amplitude equation which is of gradient 
form, and allows the explicit calculation of the functional form of the cubic nonlineari- 
ties. In particular, we find that three-wave resonant interactions contribute significantly 
to the nonlinear terms, and therefore are important for pattern selection. Minimiza- 
tion of the associated Lyapunov functional predicts a primary bifurcation to a standing 
wave pattern of square symmetry for capillary-dominated surface waves, in agreement 
with experiments. In addition, we find that patterns of hexagonal and quasi-crystalline 
symmetry can be stabilized in certain mixed capillary-gravity waves, even in this case of 
single frequency forcing. Quasi-crystalline patterns are predicted in a region of parame- 
ters readily accessible experimentally. 



1. Introduction 

The generation of standing waves on the free surface of a fluid layer that is oscillated 



vertically is known since the work of Faraday (1831) . Recently, there has been renewed 



experimental interest in Faraday waves as an example of a pattern-forming system. Rea- 
sons include the ease of experimentation due to short characteristic time scales (of the 
order of 10 -2 seconds), and the ability to reach very large aspect ratios (the ratio of 
lateral size of the system to the characteristic wavelength of the pattern) of the order of 
10 2 . By varying the form of the driving force and by using fluids of different viscosities, a 
number of interesting phenomena have been observed including the emergence of stand- 
ing wave patterns of different symmetries near onset (Christiansen, Alstr0 m & Levinsen 
1992; Fauve et al. 1992; Edwards & Fauve 1993, 1994; Miiller 1993), secondary insta- 
bilities of these patterns when the amplitude of the periodic driving force is increased 
(Ezerskii, Korotin & Rabinovich 1985; Daudet et al. 1995), and spatiotemporal chaotic 
states at even larger amplitudes of the driving force (Tufillaro, Ramshankar & Gollub 
1989; Gollub & Ramshankar 1991; Bosch & van de Water 1993; Bridger et al. 1993; 
Kudrolli & Gollub 1996). 

A numerical linear stability analysis of the Faraday wave problem has been carried out 
by Kumar fc Tuckerman (1994)| for a laterally infinite fluid layer of arbitrary viscosity. 
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The predicted values of the acceleration threshold and the wavelength at onset are in 
good agreement with experiments in large aspect ratio systems (see also Bechhoefer et 
al. 1995). In the simpler case of an ideal fluid, a classical linear stability analysis leads to 
the Mathieu equation for the interface displacement (Benjamin & Ursell 1954). On the 
other hand, the nonlinear evolution above onset, including pattern selection, secondary 
instabilities and the transition to spatiotemporal chaos are not very well understood 
theoretically. In this paper, we present a weakly nonlinear analysis of Faraday waves 
driven by a sinusoidal force in the limit of weak dissipation and for a laterally infinite 
system of infinite depth (or unbounded free surface waves). Our main focus is on pattern 
selection near onset, and our results are compared to experiments involving fluids of low 
viscosity, in containers of large depth and aspect ratio. 

Many studies of free surface waves in incompressible fluids have focused on the inviscid 
limit (see, for example, Yuen & Lake 1982; Craik 1985). Viscous dissipation, however, is 
essential for Faraday waves because it not only sets a threshold value of the driving force, 
but also affects nonlinear saturation of the parametric instability. For unbounded free 
surface waves in the limit of weak dissipation, the flow remains potential except in a very 
thin layer at the free surface (Lamb 1932; Landau & Lifshitz 1959). A common procedure 
in this limit involves the introduction of the so-called quasi-potential approximation 
(QPA) which perturbatively incorporates weak viscous effects by introducing modified 
boundary conditions for the otherwise potential bulk flow. By performing a formal 
expansion in the small thickness of the viscous boundary layer, Lundgren & Mansour 
(1988) derived a set of quasi-potential equations (QPE's) that included nonlinear viscous 
contributions for the free surface flow of an axially symmetric liquid drop. Ruvinsky, 
Feldstein & Freidman (1991) later derived a set of QPE's that contain only linear damping 
terms for two-dimensional free surface waves. In §||, we present QPE's for parametric 
surface waves including only linear viscous terms. These equations are a direct extension 
of those of Ruvinsky, Feldstein & Freidman for the two-dimensional case. They can also 
be derived in a formal expansion similar to that of Lundgren & Mansour (1988) and 
neglecting nonlinear viscous terms. Such formal derivation and the expressions for the 
nonlinear viscous terms can be found in Zhang (1994). 

Since we are interested in the surface displacement, but not in the flow field in the fluid 
interior, additional simplification is achieved by writing the three-dimensional QPE's in 
a two-dimensional nonlocal form, that involves flow variables at the free surface alone. 
Such simplification is possible because the velocity potential satisfies Laplace's equa- 
tion in the fluid interior, and is determined uniquely when its values on the boundary 
(the free surface) are known. An outline of the derivation of the weakly nonlinear two- 
dimensional form of the QPE's is given in § |2.2| by using the so-called Dirichlet-Neumann 
operator (Craig 1989). The projection of the fluid equations onto equations involving 



the coordinat es of the free surface follows the works of Miles (1977) , Milder (1977) , and 



Craig (1989) for inviscid gravity waves. These two-dimensional nonlocal quasi-potential 
equations are our starting point for our weakly nonlinear analysis of parametric surface 
waves which is presented in Standing wave amplitude equations are derived by using 
a multiple scale perturbation expansion. Pattern selection in Faraday waves near onset 



is discussed in § 3.3 



Although the relation between our results and earlier theoretical work on Faraday 
waves is discussed in more detail below, we note here two important studies addressing 
the weakly nonlinear regime above onset by Milner (1991), and Miles (1993, 1994). They 
obtained amplitude equations for inviscid flow, and then introduced weak viscous effects 
by adding damping terms directly to the amplitude equations from an energy balance 
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consideration. The resulting amplitude equations differ from ours in several respects, 
hence our results differ qualitatively from theirs, as discussed in §[| 



2. Quasi-potential equations for free surface waves 

We consider a reference state in which a quiescent and incompressible Newtonian liquid 
of density p and kinematic viscosity v occupies the z < half space, and a gas of negligible 
density occupies the z > half space. The gas phase is also assumed to have a uniform 
and constant pressure field pq. Under these conditions, the velocity distribution of the 
gas phase may be allowed to remain unknown and we can write the governing equations 
for the liquid phase only. Such a liquid-gas interface is usually called the free surface of 
a liquid (Batchelor 1967). The governing equations for the velocity field v(x,y,z,t) are 
for — oo < z < h(x,y,t), where z = h(x,y,t) is the instantaneous location of the free 
surface, 

V-v = 0, (2.1) 

<9,v+ (vV)v = --X7p + v\7 2 v + g{t)z, (2.2) 
P 

where z is the unit vector in the positive z direction, and g(t) = —go — g z (t) with go 
a constant gravitational acceleration, and g z (t) the effective acceleration caused by the 
vertical oscillation of the fluid in the Faraday experiments. The boundary conditions at 
the free surface are, 



d t h/y/l + {Vh) 2 =vn, (2.3) 

a ■ t ■ n = 2pva ■ e • h = 0, (2-4) 

b-T-h = 2pisb-e-h = 0, (2.5) 

p — fi ■ t ■ n = p — 2pvh ■ e • n = po + Tk, (2-6) 

with boundary condition, 

v — 0, as z — > — oo, (2.7) 

where a(x,y,t) and b(x,y,t) are two tangential unit vectors on the free surface, and 
h(x, y, t) = (—h x , —h y , 1)/ \Jl + (V/i) 2 is the unit normal vector of the free surface point- 
ing away from the liquid. T is the surface tension, k is the mean curvature of the free 
surface, which is given by k = V • h, and r = 2pve is the viscous stress tensor, where e 
is the rate of strain tensor with Cartesian components (i,j — x, y, z) 

The equations governing fluid motion can be simplified in the limit of weak viscous 
dissipation. In this thin viscous boundary layer, also known as the vortical 

layer, occurs near the free surface as a result of the nonzero irrotational shear stresses. 
This small irrotational shear stress drags a thin viscous layer of rotational fluid along, 
causing a small modification in the velocity field which is required in order to satisfy 



the zero-shear-stress boundary conditions (Eqs. (2.4) and (2.5)). Since the free surface 
vortical layer is thin, the basic idea of the quasi-potential approximation is to consider 
pure potential flow in the bulk that satisfies effective boundary conditions on the moving 
surface to account for weak viscous effects. 
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2.1. Three-dimensional form of the quasi-potential equations 

Let <p(x, y, z, t) be the velocity potential for the bulk potential flow. As a direct extension 
of the QPE's of Ruvinsky, Feldstein & Freidman (1991) for two-dimensional surface 
waves, the governing equations for unbounded three-dimensional surface waves read, 

V 2 (/)(a;,y,z,t) = 0, for z < h(x,y,t), (2.9) 

with boundary conditions at the free surface z — h(x,y,t) 

dth + V<j>'Vh = dz<f> + W(x,y,t), (2.10) 

dt + k^^) 2 -g(t)h + 2 V ^ = T K , (2.11) 
2 az l 

d t W(x,y,t) = 2u(d xx +dyy)d z cf>, (2.12) 

and 

d z <p — * 0, as z — > — oo, (2-13) 
where W(x, y, t) is the z-component (or the linearized normal component) of the rota- 



tional part of the velocity field at the free surface. The viscous contribution in Eq. ( 2.11 ) 
is related to the normal stress of the irrotational velocity component. We note that 
only linear viscous (damping) terms are retained in the above set of equations, which 
we shall refer to as LDQPE's. Nonlinear viscous contributions can be obtained by us- 
ing a formal perturbation expansion, similar to that of Lundgren & Mansour (1988), of 



Eq. (2.1-2.7) in the small thickness of the viscous boundary layer at the free surface. 
Details of the expansion can be found in Zhang (1994). These nonlinear viscous terms 
for three-dimensional waves are, however, algebraically too complicated to be included 
in our analysis. We note that neglecting nonlinear viscous contributions in the above 
equations is an uncontrolled approximation, motivated by the small viscosity of the fluid. 
Whether the LDQPE's are a good approximation for weakly damped parametric surface 
waves near onset is one of the central issues of this paper. Based on our analytical results 
of the LDQPE's in and the comparison of these results with experiments, we conclude 
that the LDQPE's do provide a quite good description of weakly damped and weakly 
nonlinear Faraday waves. This seems to indicate that the role of nonlinear viscous terms 
in the QPE's is not significant for pattern formation in weakly damped Faraday waves 
close to onset. 

The 2-component of the rotational part of the velocity W(x,y,t) can be eliminated. 



From Eqs. fl2.10| ) and fl2.12|) , we have, 



d t W = 2v(d XX + dyy)d t h + 2v(d XX + dyy) (V0 ■ V h ~ W) . (2.14) 

We note that vSI <f> ■ V/i is a nonlinear viscous term, and thus negligible within the 
approximation. Since W is of 0(y), the term vW in the above equation is of (D(v~ 2 ) : 
and is negligible in the weakly damped limit. Therefore we have dt(W — 2vV 2 h) = 0, or 

W(x, y, t) = 2vV 2 h(x, y, t) + W - 2jA7 2 /i , (2.15) 

where ho and W are initial conditions for h and W respectively. By considering fluid 
motion starting from at rest, we can set Wo = 0. Also since we are interested in nonlinear 
pattern formation from a nearly flat surface, ho is a small quantity. Although the term 
2v~S/ 2 ho might influence the linear growth, it is certainly negligible for nonlinear, finite 



amplitude states. Thus the boundary conditions at z = h(x, y, t) (Eqs. ( 2.10 ) and ( 2.11 )) 
now read, 



d t h = 2vV 2 h - X7 ±( f> ■ X7h + d z <f>, 



(2.16) 
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d t <t> = 2vV ± 2 cj) - i (V0) 2 + g(t)h 



Tk, 
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(2.17) 



where Vj_ = d x x + d y y. 



2.2. Two-dimensional nonlocal form of the quasi-potential equations 

For weakly nonlinear surface waves, the three-dimensional quasi-potential equations can 
be further simplified by recasting them in a form that involves only the flow variables 
on the free surface. Since the velocity potential 4> satisfies Laplace's equation in the 
bulk, it is possible to rewrite the LDQPE's as integro-differential equations involving 
variables at the free surface only. We then expand the resulting equations to third order 
in the wave steepness. Such a two-dimensional nonlocal formulation has been derived 
by [Miles (1977)1 |Milder (1977)1 |Craig (1989)| , and |Craig fc Sulem (1993)| for un-forced 
inviscid gravity waves. We extend their approach in this section to parametrically forced, 
weakly damped capillary-gravity waves. 

Let x = ( x, y), and d efine the surface velocity potential $(x, t) as <I>(x, t) = (/>(x, /i(x, t), t). 
Then Eqs. ( 2.16 ) and ( 2.17 ) can be rewritten as, 



d t h = 2vV 2 h + a/1 + {Vh) 2 d n cj), 



(2.18) 



5 t $ = 2iA7 2 $ 



i(v*) 2 



^/l + (Vh) 2 d n (t> + V$ • Vh 



2[l + (V/i) 2 ] 



+g(t)h+ -V 
P 



VT+WW 



(2.19) 



where d n <t> — n ■ (Vj_ + d z z)4>. We note that except for the normal derivative d n <p, all 
other variables only depend on the two-dimensional coordinate x. Since (/(> is a harmonic 
function, the normal derivative d n <j) at the boundary is related to the value of </> at the 
boundary. One such relation is given by the Dirichlet-Neumann operator Q(h). The 
Dirichlet- Neumann operator G(h) takes boundary values for a harmonic function and 
returns its normal derivative at the boundary with a metric pre-factor (Craig 1989; 
Craig & Sulem 1993), 



0(/i)4(x,i) = v/1 + (Vh) 2 d n 4>. 



(2.20) 



We note that Q(h) is a linear operator for $(x, t) and depends on the shape of the free 
surface z = /i(x, t) nonlocally. 

An important property of G(h) is that it has a computable Taylor expansion in powers 
of the surface displacement /i(x, t) and its spatial derivatives at /i(x, t) = 0. This Taylor 
expansion of Q(h) is useful for studying the weakly nonlinear dynamics of surface waves 
since only the first a few terms in the expansion, e.g., up to third order in /i(x, t), need 
to be evaluated. The Taylor expansion of G(h) up to order 0(h 3 ) reads, 



G(h)$(x, t) = V$ - V • (W$) - V(hV<i>) + V 



1 



hV(hV<i>) + -h z \7 z <S> 



--V 2 (/i 2 X>$) - V 2 



-h 2 V(hf><I>) + -h 3 V 2 <S> 
2 3 



V 



hV{ht>{hV$)) + i/i 2 V 2 (/iP$) - i/i 3 V 2 (£>$) + i/i£>(/i 2 V 2 $) 



(2.21) 



where T> is a linear Fourier-integral operator and is defined for an arbitrary function w(x) 
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Vu(x) = / |k|u(k) exp(ik • x)dk 
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by 



where u(k) is the Fourier transform of u(x). The operator f> is also nonlocal and is 
sometimes written as T> = \J — V 2 . By substituting the expansion for the Dirichlet- 
Neumann operator Q(h) into the boundary conditions, and consistently keeping only 
terms up to the third order in h and/or $, we have 

1. 



9t/i(x, t) = 2v\7 z h + D$ - V ■ (/iW) + -V 2 (/i 2 X>$ 



-V{hV$>)+V 

r 



hV(hf)<S>) + i/i 2 V 2 $ 



«9 t $(x, i) = 2iA7 2 $ + + -V 2 /i + - f 2 

P 2 



5 (V$)2 



(2.22) 



(2.23) 



Consideration of higher order terms is certainly not a difficulty in this formulation, but 
most of the phenomena in Faraday waves near onset should be explained within a frame- 
work that includes up to third order nonlincaritics. Equations ( 2.22 ) and ( 2.23 ) are 
two-dimensional, and they are the starting point for the analytical asymptotic anal- 
ysis presented below, and for extensive numerical studies that will be reported else- 
where (a short summary of both analytical and numerical results can be found in 
Zhang fc Vinals (1996| . 

Finally, we note that the incompressibility condition (Eq. ^J) implies that the average 
level of the surface displacement h(x,t) is constant, 



/ /i(x, t)dx = constant. 
Js 



(2.24) 



It is easy to see that J s /i(x, t)dx is indeed a constant of motion for Eq. ( 2.22 ). 

Before we proceed any further, it is useful to discuss at this point the dissipation func- 
tion approach used by other authors to obtain dissipative contributions to the equation 
of motion for weakly damped waves, and compare it with the quasi-potential approx- 
imation. As we show below, both methods already differ in the linear viscous terms 
in the dynamical equations even though, by construction, they give the same (correct) 
rate of decay of the energy for linear surface waves. As a consequence, is seems to us 
that nonlinear viscous terms obtained from an energy balance on the inviscid amplitude 
equation may not be reliable. 

It is well known that the irrotational surface wave problem in an inviscid fluid can be 
written in a Hamiltonian form (Zakharov 1968; Miles 1977). The governing equations 
for irrotational surface waves can be written in this case as, 



d t h(x,t) = 
a t $(x, t) = 



SH 



(J*(x,*)' 

SH 
~ 6h(x, t) ' 



(2.25) 
(2.26) 



where /i(x, t) and ^(x, t) are the generalized coordinate and momentum respectively, and 
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the Hamiltonian H is given by 



H = dx 



MX,t) dz ((V ± </>) 2 + (d z <P)A + \g{t)h 2 + T (yi + CV^) 2 - 1 



(2.27) 

This Hamiltonian formulation (or equivalently the corresponding Lagrangian formu- 
lation) offers a natural way to incorporate the effects of viscous damping by adding a 
dissipation function. This is done by modifying the equation for the generalized momen- 
tum $(x,t), 

SH 

5 t $(x, t) = -jr—r + Q(h(x, t), $(x, t)), (2.28) 
<m(x, t) 

where Q(h, $) is a dissipative or damping force, not invariant under time reversal, and 
often of phenomenological nature. In the case of Faraday waves, the dissipation function 
Q(h, <f>) has been determined by equating the rate of energy loss in this near-Hamiltonian 
formulation, 

^ - ^ = j "/ d *Q( h > *)9th, (2.29) 
to the decay rate of the total energy for potential flow (Landau & Lifshitz f959), 

J J dxQ(h, $)d t h = -v J J dxj ( ' V 2 (V</>) 2 . (2.30) 

It is easy to show that the linear part of Q depends only on $(x, t), and is given by 

Q(h, $) = 4i/V 2 $(x, t) + nonlinear terms. (2-31) 



If we consider the linear approximation for Q, Eq. ( 2.26| ) is modified by the addition of 



a viscous damping term 4i/V 2 $(x, t) to the RHS, while Eq. (2.25) remains unchanged. 

However, in the quasi-potential approximation (QPA), viscous damping terms appear 
in both the dth equation and the equation (Eqs. ( ]2.22|) and ( 2.23] )). This difference 



has important implications for the standing wave amplitude equations to be derived in § [| 
Recall that the viscous correction term in the dth equation is related to the rotational 
component of the velocity field at the fluid surface, while the viscous correction term in 
the dt& equation is related to the normal stress of the irrotational component. 

For linear surface waves, both approaches give the same dynamical equation for a 
Fourier mode hk(t) of h(x,t), 

/Tk 3 \ - 

d u h k {t)+Avk 2 d t h k {t)+ l— + g Q k + kg z (t)\ h k (t)=0. (2.32) 

Thus, the decay rate of the total energy for linear surface waves is the same for the two 
approaches, as expected. 

Since the calculation of the dissipation function Q can be carried out for the nonlinear 
terms as well, it would appear that the dissipative function approach is a natural way 
to incorporate nonlinear viscous terms order by order. However, since the dissipation 



function approach does not give the correct linear viscous terms (Eqs. ( 2.28 ) and ( 2.31 ) 
lead to a linear viscous term in the equation for <9 t $ equal to 4^V 2 $, whereas the linear 



viscous term in Eq. ( 2.23| ) is only 2i/V 2 <!>), there does not seem to be any a priori reason 



to trust nonlinear viscous terms. 
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3. Weakly nonlinear analysis 

Standing wave patterns of square symmetry are observed near onset in Faraday ex- 
periments of weakly viscous fluids, in containers of large lateral size, and with single 
frequency sinusoidal forcing (Faraday 1831; Rayleigh 1883; Lang 1962; Ezerskii, Korotin 
& Rabinovich 1985; Tufillaro, Ramshankar & Gollub 1989; Ciliberto, Douady & Fauve 
1991; Christiansen, Alstr0 m & Levinsen 1992; Bosch & van de Water 1993; Edwards & 
Fauve 1993, 1994; Miiller 1993). We derive next a set of coupled standing wave amplitude 
equations valid near onset that can accommodate patterns of arbitrary symmetry on a 
two-dimensional surface. The standing wave amplitude equations that we will obtain are 
of gradient form, and thus minimization of the resulting Lyapunov functional determines 
the symmetry of the most stable standing wave state. Our derivation of the standing 
amplitude equations has three novel features. The first one is the different starting point 
for the asymptotic expansion. It is based on the LDQPE's described above. Second, we 
note that there are two independent small parameters in this system, namely the reduced 
dimcnsionless driving amplitude e, which is also the distance away from threshold, and 
the viscous damping parameter 7 (to be defined below). A double perturbative expansion 
for these two small parameters is necessary. Solutions of the linearized quasi-potential 
equations are obtained by performing a perturbative expansion for the small damping 
parameter or the driving amplitude / (to be defined below) . The linear solutions contain 
the primary mode of the fluid surface with a frequency half of the driving frequency 
as well as its higher harmonics. These higher harmonic terms are proportional to the 
driving amplitude / or its powers (/", with n = 2, 3 • • • .). The nonlinear interaction of 
these higher harmonic terms with the primary mode provides a novel amplitude-limiting 
effect for the parametric surface wave system. This effect is important for the nonlinear 
saturation of the surface wave amplitude in weakly dissipative systems. 

The third feature is related to three-wave resonant interactions in capillary-gravity 
surface waves. Although quadratic terms are prohibited by symmetry in the standing 
wave amplitude equations that we derive, three- wave resonance (triad resonance) plays a 
crucial role in pattern selection. Both three- and four-wave resonant interactions among 
capillary-gravity waves arc well known and well studied. The importance of three-wave 
resonant interactions to pattern selection in Faraday waves, however, has been largely 
overlooked. As we show later, the resonant interactions between two linearly unstable 
standing wave modes and a linearly stable wave mode strongly affects four- wave nonlinear 
interactions, and thus the coefficient of third order nonlinear terms in the amplitude 
equations. 



Previous theoretical work by Milner (1991) involved the derivation of a set of cou- 



pled traveling wave amplitude equations for inviscid parametric surface waves, to which 
viscous damping terms were added by an energy-balance-consideration, equivalent to 
the dissipation function approach describe above. He concluded that nonlinear viscous 
damping terms in the dissipation function play a major role in pattern selection. We 
disagree with his conclusion for four reasons: (i) As discussed earlier, the dissipation 
function approach does not give the correct linear viscous terms, so it is doubtful that it 
will introduce the correct nonlinear damping terms; (ii) linear viscous terms in the fluid 
equations can contribute to nonlinear damping terms in the amplitude equations, while 
such contribution is absent in Milner's phenomenological consideration of viscous effects; 
(iii) an amplitude-limiting effect of the driving force did not appear in Milner's analy- 
sis since he used a zeroth order linear solution for the parametric instability; and (iv) 
Milner did obtain triad resonant interactions in his calculation, but by not taking them 
into account explicitly, he overlooked their effect on pattern selection. As was recently 
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suggested by Edwards & Fauve (1994), we will show that triad resonant interactions play 



an important role in pattern formation of Faraday waves in weakly viscous fluids. 

3.1. Solutions of the Linearized Equations 

As is well known, the linearized problem of parametric surface waves can be reduced to 
the damped Mathieu equation, and the Faraday instability corresponds to the subhar- 
monic resonance of the equation. For the case of a sinusoidal driving force, the effective 



acceleration g(t) in Eqs. ( 2.22 ) and ( 2.23 ) can be written as g(t) = —go — g z sinfii, where 
go is the constant acceleration of gravity, and Q and g z are the angular frequency and 
the amplitude of the driving force respectively. We now choose 2/0 = l/wo as the unit 
of time and l/fc as the unit of length with fc defined by uj 2 = goko + -k 3 ,. We also 
choose the unit for the surface velocity potential $ to be ujo/ko- We further define a 
dimensionless linear damping coefficient 7 = 2vko/uj , Go = .go^o/^o, r = Tkof(puJo), 
and the dimensionless driving amplitude / = g z ko/ '(4uq). Note that Go + To = 1 by 
definition. 



By linearizing the quasi-potential equations and boundary conditions (Eqs. (2.22) and 
( |2.23 )) with respect to the surface displacement h, and the surface velocity potential 



and taking the Fourier transform with respect to the spatial coordinate, one obtains in 
a standard way, 

d tt h k + 2jk 2 d t h k + (G k + T k 3 + 7 2 fc 4 + 4/fcsin2t) h k = 0, (3.1) 

k$ k = d t h k + jk 2 h k . (3.2) 



Equation (3.1) is the damped Mathieu equation for h k . We now seek analytical solutions 
of the above equations perturbatively. We introduce a small parameter rj (77 <C 1) such 
that 7 = 7 o?7, and / = / ?7, where 70 and f are assumed to be of O(l). For weakly 
dissipative fluids, i.e. 7^1, and to be consistent with the quasi-potential approximation 
discussed in last section, we will neglect the term proportional to 7 2 in Eq. (3.1). At 



subharmonic resonance, we have uj 2 (k) = Gok + Tok 3 = 1, which of course implies 
k = 1 (note that Go + To = 1 by definition). We then consider an expansion for the 
wavenumber k near subharmonic resonance as, k = 1 + Akr/+ ■ ■ ■ . Above (and near) the 
onset of subharmonic resonance, we expect the amplitudes of h k to grow in time but in a 
slower time scale than that for the subharmonic oscillation. In the following we assume 
that the slow time is T = T)t, and seek solutions perturbatively as power series in 77, 



h k = h^(t,T)+ v hl"(t,T) + 



At O(?7 ), one has, 



h { " ] (i, T) = A k (T) cos t + B k (T) sin t, 
$i 0) (i,T) = -A k (T)smt + B k (T) cost, 

where A k (T) and B k (T) are arbitrary functions. At O^q 1 ), a standard solvability condi- 
tion appears, 

d T A= (/o- 70)^+^ (Go + 3r )AA;S, (3.3) 
d T B = -(/„ + 7o )B - i(G + 3r )AfcA (3.4) 
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By substituting A and B oc e CTt , one finds 



°± = ~7o ± y /o 2 - (^(Go + 3r )) 2 . (3.5) 

Exactly at subharmonic resonance (k = 1 or Afc = 0), the growing mode M + oc A and 
the decaying mode M_ oc B. In this case, the linearly growing eigenmode above onset is 
given by, 

/ x N 

'"i az i i 

' j=i 

JV 

c 

3=1 

where we have assumed that the standing waves consist of an arbitrary discrete set of 



h(x, t) — ycos t + — sin 3t + ■ ■ ■ j ^ ^Aj (t) exp (ikj ■ xj + c.c. , (3.6) 
$(x, t) = (-smt + fcost + ^-cosM H j ^ ^Aj(t) exp Tifej ■ + c.c. , (3.7) 



wavevectors in the two-dimensional space. When / = 7 (at onset), Eqs. ( J3.6| ) and (3.7) 
are the linear neutral solutions, which is the basis of a weakly nonlinear analysis for our 
problem. It is important to note at this point that we have kept in the linear solution 
terms proportional to / (or 7 since / = 7 at onset). These terms will be crucial for 
obtaining the correct cubic term in the amplitude equations, and had not been included 
in previous studies. Terms proportional to higher harmonics do not contribute to the 
standing wave amplitude equations to the order considered. 

3.2. Standing Wave Amplitude Equations 
We seek nonlinear standing wave solutions of Faraday waves near onset (e = (/ — 7) <C 



1) in this section. We expand the two-dimensional quasi-potential equations (Eqs. (2.22) 



and (2.23)) consistently in e 1 / 2 with multiple time scales, 

h(x, t, T) = e^/nfo t, T) + eh 2 + e 3/2 h 3 + ■■■, (3.8) 

$(x,i,T) = £ 1/2 $ 1 (x,i,T)+e$ 2 +e 3/2 $ 3 + ---, (3.9) 

where the slow time scale T = st, and is not related to the slow time in the previous 
section. The scaling of the amplitudes and the slow time T with e is formally determined 
by the ultimate consistency of the expansion, and in particular by the balance of terms 



in the final form of the standing wave amplitude equation (Eq. (3.34)). We define the 
following linear operator, 



C - ' Go - rtv 2 + V 4 7 sin It d t - 7V 2 ) ' i: :i ' ' ° ' 



On substituting Eqs. (3.8) and ( pj| ) into the quasi-potential equations we have at 0(e 1 ^ 2 ) 



£ ft )=»■ 

The above equation is the same as the linearized quasi-potential equations discussed in 
the last section except the driving amplitude / is replaced by the damping coefficient 7. 
Thus, the solutions for hi and $1 are just the linear solutions found in the last section. 
For simplicity, we will neglect the linearly stable mode B of the linear solutions. Then, 
h\ and $1 are the linear neutral solutions, 

N 

hi = (cost+ j sin 3^ \Aj(T) exp (ik 3 ■ xj + c.c. , (3.11) 

.7=1 



Pattern formation in weakly damped Faraday waves 

N 



$1 = 

At 0(e), we have 



3 7 

sin t + 7 cos t H cos 6t 

4 



[Aj(T)exp(ikj x 



11 
(3.12) 



3=1 



4> 2 



(3.13) 



-V ■ (%V$i) -X> 

l(P$ 1 )2- 1(V$!) 2 

The above equation represents two coupled equations for h 2 and $2 . It is easy to obtain 
an independent equation for h 2 from Eq. (3.13), which reads, 

d tt h 2 - 2jV 2 d t h 2 + (Go - T V 2 )t>h 2 + 4 7 /i 2 sin2t 
V • (/iiV$i)+2> (hiD$x) 



= X -V [(P$i) 2 - (V$r) 2 ] + (7V 2 - d t 

3,1=1 { 
- 7 sin 2t 



2(l + c j/ )-cos2t 



1 + Cj 



(5/2 + c,-,) ( 1 + c 3l - ^2{l + Cjl )\ + -±£zi^2(l + CjJ ) 



A, A; exp 



fc/) ■ X + c.c. 



N 

+ E< 

-7 sin 2£ 



1 — Cj 



2(1 - Cji) - cos2i 



(5/2 - ty) (l - Cjl - ^2(1 - Cj ,)) + ^-^^2(1 -c jt ) 



AjA* exp I i(fcj — fcj) • x) + c.c 



(3.14) 



where Cj-j = cos 9 ji = kj ■ ki and we have neglected terms that are proportional to 7 2 . 
We note that there are no terms on the RHS proportional to cost or sint, that would 
introduce a secular variation in the solution. Therefore, there is no solvability condition 
for the amplitudes Aj at this order. There are, however, resonant interactions due to 
certain terms on the RHS. 

The particular solution h 2 of Eq. (3.14) can be written as 

N 



3,1=1 



AjAi exp I i(kj + ki) • x) + c.c 



+Hj^-i(t) AjAi exp [i(kj — ki) • x) + c.c 



(3.15) 



where Hji(t) is an unknown function to be determined, and Hj-i is defined by replacing 
Cji with Cj y -i in Hji. On substituting the above form for h 2 into Eq. ( |3.14| ), we have the 
following equation for Hji (t) , 



d tt H 3l + 27^/2(1 + crfdtHfl + [Go + 2r (l + c j i)] y /2(l + c i i)H j i 



F ( ,P cos 2t + Ff, 2) sin 2t 



(3.16) 



where F^/ and F^f' are proportional to AjAi. Only terms that are relevant to the 
resonance are written out on the RHS of Eq. ( 3.16j ). Equation ( 3.16| ) looks very much 
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like the equation for an additively forced harmonic oscillator with friction. When the 

"natural" frequency of the "oscillator", [(Go + 2F (1 + c j i))- s /2(l + c jZ )] 1/2 , equals the 
driving frequency, resonance occurs. This condition reads, 

[G + 2r (l + c j ,)]^/2(l + c J -,) -4. (3.17) 



Due to the nonzero damping coefficient, 2jy/2(l + Cji), this resonance results in a finite 
value for Hji that is inversely proportional to the damping coefficient. We note that the 
parametric forcing term hi sin 2t in Eq. ( |3.14| ) is not directly relevant to the resonant 
interaction. 



The values of 9ji (cji = cos9ji) that satisfy the resonance co nditio n (Eq. (3.17)) as a 
function of r are shown in Fig. 0(a). Since the RHS of Eq. ( 3. 16|) is proportional to 



AjAi, there are three waves involved in this resonance, namely, standing wave modes A 
and Ai, and mode B with wavevector kj + ki, as shown in Fig. |l|(b). Therefore, Eq. (3.16) 



describes a three-wave resonant interaction. Note that the wavenumber for mode B is 
away from the critical wavenumber ko = 1, and thus mode B is a linearly stable mode. 
For To < 1/3, triad resonance is not possible. For Tq — 1/3, wavevectors of the three 
resonating waves are in the same direction (6^ =0). As T is further increased, 6^ 

also increases. For purely capillary waves (Tq = 1), 9^ reaches the maximum value of 

Cjl = 2 1 / 3 - 1 or 9^ « 74.9°. 

Resonant interactions among surface capillary-gravity waves in general have been stud- 



ied since the pioneering work by Wilton (1915) (for a recent review see Hammack & 



Henderson 1993). Wilton found that a Stokes expansion in powers of the wave ampli- 
tude became singular at a wavenumber k — y // gop/2T for inviscid capillary-gravity waves 
in two spatial dimensions, which corresponds to the triad resonance discussed above for 

r = i/3. 

This special case of resonant interaction is often termed second-harmonic resonance 
and the corre sponding capillary-g ravi ty waves are called Wilton 's ripples. Experimen- 



tal studies by McGoldrick (1970) , and Banerjee fc Korpel (1982) on the triad resonance 



in capillary-gravity waves have verified the function 9^ (see Fig. |l|(a)) quantitatively 
(Hogan 1984). 

The relevance of triad resonant interactions at second order to pattern formation in 
Faraday waves has been largely overlooked by previous studies. We want to emphasize 
here that the effect of triad resonant interactions on pattern formation in parametric 
surface waves can be intuitively understood. Let us consider a situation in which two 
linearly unstable standing waves Aj and Ai with their wavevectors separated exactly the 

(r) 

resonating angle 9 -^ grow from small amplitudes in the linear regime. When they enter 

the nonlinear regime, a mode B with wavevector kj + ki is created as a results of quadratic 
nonlinear interaction of Aj and Ai . The amplitude of B will become very large because 
of the resonance and weak damping. For inviscid fluids, the amplitude of B will increase 
without limit. Since the parametric forcing pumps energy into the surface wave system 
through unstable modes Aj and Ai at a rate determined by the supercriticality e, the 
energy required for the growth of mode B will come from reductions of the amplitudes in 
Aj and Ai, but not directly from the external parametric force. As a result, the modes 
Aj and Ai will have smaller amplitudes than other unstable modes that do not satisfy 
the triad resonance condition. In other words, small perturbations with components 
close to the critical wavenumber will grow exponentially. When the amplitudes of these 
modes become large enough, a nonlinear selection process takes place such that modes 
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(r) 

with their wavevectors separated by the resonating angle 9^ are less favored or avoided. 
The above argument is relevant for both inviscid or weakly viscous fluids. For fluids of 
high viscosity, the influence of triad resonant interactions becomes smaller because of the 
large damping. 

Since there is no solvability condition at second order, there will be no quadratic 
terms in the amplitude equations to be derived. Thus, the above described influence 
of triad resonant interactions must appear through cubic nonlinear terms, which repre- 
sents four- wave resonant interactions among four linearly unstable standing wave modes. 
Specifically, in particular the coefficient of the cubic nonlinear terms in the amplitude 

(r) 

equations are expected to have a large positive peak at the resonating angle 9 ^ as we 
show later in this section. When viscous damping effects are neglected, the peak shifts 
toward +00. Such divergence of the nonlinear interaction coefficient was encountered 



by Milner (1991) in his analysis of weakly damped parametric surface waves because he 
neglected the contribution to the amplitude equation from the second order solution. As 
a consequence, he failed to realize the relevance of such a divergence to pattern se lec- 



had 



tion. In contrast, and in agreement with our calculations, Edwards fc Fauve (1994) 
suggested recently that triad resonance could be important for pattern selection. 

As we show later, the occurrence of standing wave pattern of square symmetry in cap- 

(r) 

illary Faraday waves is closely related to triad resonance withfl}^ = 74.9°. By increasing 

the gravity wave component in capillary-gravity waves, the resonating angle 9 ^ becomes 
smaller and the nonlinear interaction coefficients at cubic order changes accordingly. For 
To ~ 1/3, 6jp is close to zero and we will show that standing wave patterns of square 
symmetry become unstable and hexagonal, triangular, or quasicrystalline patterns can 
be stabilized. 



The solution of the corresponding homogeneous equation of Eq. (3.14) will have the 
same form as the solution at order 0(e 1 ^ 2 ), and thus can be absorbed into the solu- 
tion at 0(e 1 ^ 2 ). As a result, we are only interested in the particular solutions to the 
inhomogeneous equation for /12. The particular solutions are, 



N 

E 

3,1=1 



(aji + [3ji cos 2t + jSji sin 2t) A 3 Ai exp (i{k 3 + ki) ■ xj + c.c. 
(<x,7 -I- Pji cos2i + jSji sin2i) A 3 A\ exp (i(kj — fc;) • x^j + c.c. 1, (3.18) 



where 



and 



31 4[G + 2r (l + Ci/ )] G + 2r (l + c, i )' [ ! 
- (l + Cji - ^y/2(l + Cjl )) (D M - 8 7 2 M j7 ) + 8 7 2 (1 + C]l )N 3l 

Pjl ~ 64 7 2(1 + CjJ ) 2 +^-8 7 2^M j7 ' (3 - 20) 
8(1 + Cjl ) (l + Cjl -^ ^/2(l + c 3l )) + D 3l N 3l 

6jl ~ 64 7 2 (l + c, i ) 2 + ^- 87 2 ^M,/ ' (3 ' 21) 

= G + 2r„(i + c J ,)' (3 - 22) 

Dji = [Go + 2T (l + cji)] ^2(1 + c jt )- 4,(3.23) 



14 



■Cji 



W. Zhang and J. Vinals 

1 + Cfl 



1 + Cjl - yj2(l + Cj i)j + — g^y / 2(l + cj0 + (1 + c,-i)AOj (3.24) 



cijj, and can be obtained by replacing Cji with —Cji in the expressions for a 
j3ji, and Sji respectively. The triad resonance occurs when Dji = 0. We have retained 
terms that are proportional to 7 2 in the expressions for ctji, (3ji, and Sji since these terms 
become important when Dji is small or zero, i.e., at triad resonance. The factor 7 2 in 
these terms will either cancel to give terms of 0(1), or partially cancel to give terms of 
0(1/7). Although we still keep them when away from resonance, these terms are very 
small for weak damping, and thus should not affect the consistency of the perturbation 
expansion. 



From Eq. (3.13), we have 

P$ 2 = d t h 2 - 7 V 2 /i 2 



(3.25) 



On substituting the expressions for hi, and h 2 into the above equation, we obtain 
the following expression for $2, 



"l (jUji + jVji cos 2i + Wji sin 2t) AjAi exp (i(kj + fcj) • xj + c.c. 
+ (jUji + jVji cos 2t + Wji sin 2t) AjAi ex P ~ fa) ' x ) + c - c - r i (3-26) 



JV 

E 



where 



+ an 



3 
4 



+ ft/ - 



Wji = 



2(1 + Cjl ) 
7^2(1 + ^) 



^(l + c,-,), 

y/2(l + Cji y 



(3.27) 
(3.28) 
(3.29) 



V2(l + Cj 

and , Wjj, and can be obtained by replacing Cji with —Cji in the expressions for 
Uji, Vji, and w^; respectively. We also note that Vji and Wji are not singular at Cji = — 1 
since the factor in deno minat or y /2(l + c,;) will be canceled by the same factor in the 
numerator in Eqs. ( 3.28| ) and ( |3.29] ). Similarly, Vji and Wji are also not singular at Cji = 1. 

The amplitude Aj's are not determined yet at 0(e), and thus it is necessary to continue 
the expansion to higher orders. At 0(e 3 / 2 ), we have 



C 



h 3 



-d T hi - V • (/iiV$ 2 + h 2 V$i) - V (h{D<S> 2 + h 2 T>$i 
-dr$i - 47/11 sin 22 + X>$iX>$ 2 - V$i • V$ 2 




/iiV 2 $i 



V 2 l/l 2 ^*! 



tV- (V/ii(Vft 2 )) 



(3.30) 



We have found that it is convenient to obtain the solvability condition at this order from 
the independent equation for h%, which can be easily obtained from Eq. ( 3.30| ), and it 
reads, 

d tt h 3 - 2jW 2 d t h 3 + (G -r V 2 )m 3 + ijsin2Wh 3 = -&r ((7+&)>H+*i) 



-47/11 sin 2t + (7V 2 - d t 
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V 



+V I -h}$i - hiD(hi$i) \ - -V 2 (/i^i 
$iP$ 2 -V$i -V$ 2 + §\hx - §ib{h x $x) - — V- (V/ii(V/ii) 2 ) 
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(3.31) 



The Fredholm alternative theorem requires that the RHS of Eq. ( 3.31 ) be orthogonal 
to any of the independent solution of its adjoint homogeneous equation. The solvability 
condition reads, 



dt 



d@ j RHS 3 h j (x,t) = Q, 



(3.32) 



where RHS3 stands for the RHS of Eq. (3.31), and hj(?t,t) is either one of the following 
two independent solutions of the adjoint homogeneous equation, 



(cos t + — sin 3t 
V 4 



exp 



(±ikj ■ x) 



sint 



■ cos 3t 



exp (±ikj ■ x^ 



where kj is a unit vector in arbitrary direction, and Oj = kj ■ x. Since both RHS3 and 



/ij(x, t) have a temporal period 2n, we have taken the interval of integration in Eq. (3.32) 



to be 2tt. Also because of the periodicity of RHS3 and hj(x.,t), it is only necessary to 
consider the term proportional to cost or sint in /ij(x, t). It turns out that the relevant 
solvability condition is obtained from the second solution of hj(x,t). If we had also 
considered the linearly stable mode Bj in this analysis, another solvability condition 
would be obtained from the first solution of h(t). In summary, the solvability condition 
for Eq. (3.31) is that the coefficient of sintexp (±ikj ■ x") term in RHS3 equals zero. 



In what follows, we will collect terms from the RHS of Eq. ( |3 .31 ) that are relevant to 
the solvability condition in a tedious but straightforward calculation. These terms can 
be written as 



RHS = 2smtY J < 

3=1 



/ /28 + 9r n 3 1 > . . l2 
d T Aj - 1 A j + 7 I 64 + 2a n + -fa - -Sjj ) \A 3 1 2 



A' 



J2 9(cji)\Ai\ 2 U 

l=l(l=£j) ) 



exp 



(ikj ■ x) 



c.c. } + 



where 



<?M = ^ (1 + 24) + I U - yj2(l + Cjl ) - y/2(l- Cjl ) 



+ ( 1 + Cj i - \l 2(1 + Cj i) ) ( -wji - Vji 



+ I 1 - c jt - y 2(1 - cjt) j \ -Wji - v,i 
+ < I <■„ ) ( 2n ;/ + 1^., - -«,,) + (1 - cji) \2a 3l + - ^7 



(3.33) 
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The solvability condition therefore reads, 
dA 



dT ' 3 



N 

75 (l)L4,| 2 +7 Yu dicflW 



A h (3.34) 



where j = 1, 2, • • • , N and 



5 (1) = — ^ + 2a n + -(3 n - -6 n . (3.35) 



+ loLii H — da < 

64 JJ 8 Hn 2 

Equation ( |3.34 ) is the coupled set of standing wave amplitude equations (SWAE) for 
Aj, which is the central result of this weakly nonlinear analysis. The generic form of 
the above set of amplitude equations is of course quite general and has been derived 
for a number of different physical systems (Cross & Hohenberg 1993; Newell, Passot 
& Lega 1993). The behavior peculiar to each system stems from the functional form 
of the nonlinear interaction coefficients g(l) and g{cji), and from the time constant To 
(t = I/7, in the case of Faraday waves.) 

Before we look at the quantitative details of the nonlinear coefficients, g(l) and g(cji), 
we have the following comments on the SWAE's. 

(i) The exclusion of quadratic nonlinear terms in the standing wave amplitude equa- 
tions, which is related to the absence of solvability conditions at 0(e), is a consequence of 
the requirement of sign invariance of the SWAE's (Aj — > — Aj). Subharmonic response 
of the fluid surface to the driving force f sm(2uJot) implies /i(x, t + tt/ujq) = —h(x.,t), 
where h is a linear unstable mode given by Eq. (^J). We note that a sign change of 
the amplitude Aj is equivalent to a time displacement in a period of the driving force, 
t —* t + tt/loq. Because of the invariance of the original fluid equations under such a time 
displacement, the amplitude equation of Aj must be invariant under a sign change in Aj. 

(ii) The coefficients of cubic nonlinear terms, 7 <?(1) and 7 g(cj;), are proportional to 
the linear damping coefficient 7. This result is qualitatively different from that obtained 



by |Milncr (1991)| . He also d erived a coupled set of standing wave amplitude equations of 
the same form as Eq. ( |3.34 ) . Although his nonlinear coefficients are also proportional to 



7 (the dissipation function Q is after all proportional to the kinematic viscosity v), they 
result entirely from nonlinear viscous terms in the dissipation function. In fact, linear 
viscous terms did not contribute at all to third order. The appearance of nonlinear 
terms proportional to the linear damping coefficient 7 in the SWAE's in our approach 
is, however, no surprise. In general, a parameter that appears in the coefficients of the 
linear terms in the original equations can appear in the coefficients of nonlinear terms 
of the amplitude equations for that system. For example, the nonlinear interaction 
coefficient in the amplitude equation for Rayleigh-Benard convection is a function of 
the Prandtl number Pr, although Pr appears only in the coefficients for linear terms 
in the Boussinesq equations (see e.g., Cross 1980). Throughout this paper, nonlinear 
terms that are proportional to the kinematic viscosity v in the fluid equations (or the 
quasi-potential equations, Eqs. ( [2.9] )-( |2~i~5| )) are termed nonlinear viscous terms, and 
the nonlinear terms proportional to v or 7 in amplitude equations for Faraday waves 
are termed nonlinear damping terms in order to avoid confusion due to terminology. 
The validity of the quasi-potential equations with only linear viscous terms relies on 
the assumption that nonlinear viscous terms do not have significant effect on pattern 
formation in Faraday waves. A check of the validity can only be provided by comparing 
our results to experimental studies. 



(iii) There are contributions to the nonlinear terms of Eq. (3.34) from the parametric 



driving force. These contributions are proportional to the driving amplitude /, but they 
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appear in Eq. (3.34) together with the contributions from the linear viscous terms in 



the quasi-potential equations since we have set / — 7 in the linear solutions (Eqs. (3.11) 



and (3.12)). The contributions from the parametric driving force are directly related to 
the higher harmonic terms proportional to / in the linear solutions for h and <f>. This 
contribution provides the amplitude-limiting effect by the driving force. This amplitude- 
limiting effect results from the nonlinear interaction of these higher harmonic terms 
with the primary mode, which has half of the driving frequency. Such nonlinear in- 
teractions produces terms that are out of phase by tt/2 with the primary mode, and 
thus possibly damps or limits the wave amplitude. An important point is that this 
amplitude-limiting effect results from non-dissipative terms in the governing equation. 
Therefore, such an effect is also important for inviscid systems. This amplitude-limiting 
effect by the driving force, to our knowledge, has not been identified before. A yet 
different amplitude-limiting effect als o through non-dissipative terms was studied by 



Zakharov, Lvov fc Starobinets (1975) in the context of parametric spin-wave systems. 
Zakharov et al. studied parametric spin-wave instabilities, and considered the deviation 
of the phase of the excited spin waves from the optimum phase as the major nonlinear 
mechanism which limits the parametric instability. However, we agree with Cross & Ho- 
henberg (1993) that this "de-phasing" effect of the parametric mode is small (of higher 
order) close to threshold. 

Due to the mode interference occurring exactly at Cji = 1 (J = I or 9ji = 0), g{cji) is 
not a smooth function of Cji . It is easy to show that 

.9(1) = ±g( Cjl -> 1). (3.36) 

We also note that 

9(cji) = g{-Cji). (3.37) 
This symmetry is an obvious requirement for standing wave amplitude equations since 
it is equivalent to have two standing waves separated by angle 9 or by n — 9. 

An additional issue concerns the transformation of the equations of motion under time 



reversal, and the related question of the existence of cubic terms in Eq. ( 3.34 ) in the limit 
of a Hamiltonian system (i.e., in the absence of dissipative contributions to the equations 
of motion). Under time reversal, the variables of interest transform according to t — > —t, 



kj — > — kj and Aj(t) — > A*(—t). We also note that the Hamiltonian Eq. ( 2.27 ) depends 
explicitly on time, and that given our choice of driving force proportional to sin2i, 
the trajectory of the system under time reversal is invariant only if, in addition to the 
transformation rules given above, / — > — /. A similar situation arises in systems with 
applied magnetic fields, or for Coriolis forces, in which either the magnetic field or the 
angular velocity must be reversed. In the absence of viscous dissipation (7 = 0), the only 
nonlinear terms in the SWAE's come from the driving force, and are proportional to the 
driving amplitude /. These terms can be written as, 



-/ 



N 

<?(l)L4,f + J2 9{c 3l )\M 2 

1 = 1(1^3) 



A,. 



Thus, nonlinear terms proportional to / do change sign under time reversal, and are 
allowed in the equations of motion for the amplitudes Aj . This conclusion is not trivially 
related to having chosen a sine function as the driving force. Had we chosen / cos(2i) as 
the forcing term, the linearly growing modes would be a combination of Aj and Bj , the 
algebra for the derivation of the SWAE's becomes more tedious, but the final conclusion 
remains the same even though, in this latter case, / would be formally chosen to be 
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invariant under time reversal. On the other hand, for an autonomous system in which the 
Hamiltonian does not explicitly depends on time, if the cubic term in Eq. ( |3.34 ) is entirely 



of Hamiltonian character, then invariance under time reversal implies g(cji)\Ai \ A 
0, for any arbitrary set of amplitudes. This equality is satisfied if g(cji) = —g(—Cji), 
where Cji involves the interaction between modes kj and ki, and — Cji between kj and 
— ki (see Fig. |j]and Cross & Hohenberg (1993)). This symmetry together with Eq. ( 3.37 ) 
would imply that g(cji) = in this case. Therefore, saturation of the wave amplitude 
would have to occur cither through weak nonlinear dissipative effects, or higher order 
terms (e.g., nonlinear frequency detuning terms). 

Figure || shows the self-interaction coefficient g(l) as a function of Tq for two different 
values of 7 = 0.1 and 0.02. The maxima in the curves around Tq = 1/3 are the result 

(r) 

of the triad resonant interaction of standing waves in the same direction: 9^ — (the 
second-harmonic resonance), that affect g(l). Since the resonance is less damped for 
smaller values of 7, the value of the coefficient g(l) varies inversely with damping. 

Since g(l) > 0, we can rescale the amplitude as Aj = J g(l)Aj. We have the following 
standing wave amplitude equation for the scaled amplitude, 



IdAj 
7 dT 



N 

;=i(¥j) 



A h (3.38) 



where g{cji) = g(cji)/g(l). From Eq. ( 3.36 ) and symmetry for ±0,7, we also have g(cj 
±1)=2. 



The nonlinear interaction coefficient g(cji) in Eq. ( 3.38| ) (we have suppressed the tilde 



since we will only refer to the scaled nonlinear coefficient in what follows) depends on 
the dimensionless fluid parameters r (or Go) and 7. Figure |3| shows the function g(cji) 
for four different values of the damping coefficient 7. The maxima in g(cji) around 
Cji = 0.26 (0ji — 74.9°) correspond to triad resonance for purely capillary waves (see 
Fig. |l|). Even for relatively large values of the damping coefficient 7 = 0.2, the influence 
of the triad resonance on the g{cji) curve can still be seen, but becomes much weaker. 
An important feature common to all the curves in Fig. [| is that there is a minimum of 
g(cji) at Cji — 0, and g(0) < 1. It is also interesting to compare the differences in g(cji) 
for purely capillary waves (ro = 1), purely gravity waves (ro = 0), and mixed gravity- 
capillary waves (0 < To < 1) since the triad resonant interaction strongly depends on 
the value of To . Figure ^ shows the function g(cji) for four different values of damping 
coefficients (the same as in Fig. ||) for purely gravity waves (ro = 0). We note that 
because of the absence of triad resonant interaction, the variations among the g(cji) 
curves for different values of 7 are quite small. The curves still have minima at cji = 0, 
but the minima are much flatter. 

Since the triad resonant interaction occurs among waves with their wavevectors in the 
same direction when To = 1/3, we finally examine the function g(cji) for this case. Figure 
U shows g(cji) for four different values of the damping coefficient (the same as in Fig. 
for capillary-gravity waves with T = 1/3. We observe that g(cji) is very flat and reaches 
very small positive values for small values of 7. These facts are a consequence of the 
second-harmonic resonance, since the value of g(l) becomes very large (see Figure |]) for 
this special case of triad resonant interaction. For relatively large values of the damping 
parameter, e.g., 7 = 0.20, the effect of the second-harmonic resonance is much smaller. 
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3.3. Pattern Selection Near Onset 

Equation (3. 38) is of gradient form 1/jdr Aj = —dT/dA*, with Lyapunov function T 
given by, 

N N I N 

3=1 J'=l \ i=l(¥j) 



Since 



the only possible limiting cases of such a dissipative system, in the limit T — > oo, are 
stationary states for the amplitudes Aj . Only the states which correspond to local minima 
of the Lyapunov function are linearly stable. 



Apart from the trivial solution of Aj = for j = 1, ■ ■ ■ , N, Eq. (3.38) has a family of 
stationary solutions differing in the total number of standing waves N for which Aj ^ 0. 
By considering the case in which the magnitudes of all standing waves are the same, i.e. 
\Aj \ = \ A\, Eq. ( 3.38) ) has the following solutions, 



JV 



-1/2 



1^-1 = 14= 1+ £ 9{cji)\ . (3.41) 
The values of Lyapunov function for these solutions are, 

1 l + 2Zi=i(i^j)9{Cji) 

We note that the greater the square of the amplitude, the lower value of the Lyapunov 
function. Also the larger the values of g(cji) for a standing wave pattern, the larger 
value of the Lyapunov function. In particular, if an angle separating the wavevectors 
of two standing wave modes of the pattern satisfies the triad resonant condition, the 
corresponding g(cji), and thus the value of the Lyapunov function, will be large. There- 
fore, such patterns will not likely to appear. This result is consistent with our intuitive 
understanding of the role of the triad resonant interaction, i.e., the system tries to avoid 
pairs of standing waves with their wavevectors separated by an angle satisfying the triad 
resonant condition. 

For N — 1 (parallel roll solution), T\ = -\. For N = 2, we have either square 
(ci2 = 0) or rhombic (ci 2 7^ 0) patterns with T2 = —1/(1 +.9(212))- By considering 
regular patterns^ only, for N = 3, we have either hexagonal or triangular patterns, which 
have the same value of the Lyapunov function 

T 3/2 



1 + 5(1/2) + .g(-l/2)' 



We first consider square patterns for TV = 2. If g(0) < 1, we have T<i < T\ — — A. As 
shown in Figure [| we indeed have g(0) < 1 for the interesting parameter range of r 
and 7. Therefore standing wave patterns of square symmetry always have lower values 
of Lyapunov function than parallel roll patterns for weakly damped parametric surface 
waves near onset. 



f By regular patterns, we mean pattern structures for which the angle between any two 
adjacent wavevectors kj and kj+i is the same and amounts to n/N. 
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A32 = Tz — T2 — 



In order to compare the values of Lyapunov function for square patterns with that of 
hexagonal or triangular patterns, we compute the value of 

1 + g(l/2) + g(-l/2) - |(1 + g(0)) 
(l + .g(0))(l + fl (l/2) + .g(-l/2)) ' 

which is plotted in Figure ^. 

For 7 = 0.1, we have A32 = J-~3 — > for all values of To, and thus standing 
wave patterns of square symmetry also have lower values of the Lyapunov function than 
hexagonal/triangular patterns. We also note that the difference between Tz and T2 
becomes smaller for smaller values of Tq. For 7 = 0.02, we still have A32 = T?, — T2 > 
for capillary waves, but near the second-harmonic resonance (ro = 1/3), we have T$ < 
T2, i.e., hexagonal/triangular patterns have lower values of the Lyapunov function than 
square patterns. 

Regular patterns for N > 4 are two-dimensional quasicrystalline patterns (or quasi- 
patterns (Edwards & Fauve 1993, 1994)). A quasipattern has long-range orientational 
order but no spatial periodicity, thus analogous to a quasicrystal in solid state physics 
( Bhechtman et al. (1984) ). For N — 4, the value of the Lyapunov function for an eight- 
fold quasipattern is, 

T = ? 

l + g(V2/2)+g(0)+g(-V2/2Y 

We are interested in any parameter range in which T± has lower value of the Lyapunov 
function than Ti and The most possible parameter range is certainly near the 
second-harmonic triad resonant interaction with very weak damping. We thus compute 
the values of 



A43 = T± — T% — 

and 



(1 + g(0)) + 3g(V2/2) - 2(1 + 20(1/2)) 



((l + 2 5 (l/2))(l + .g(0) + 2.g(V2/2)) 
2 fl (V2/2)- (1+0(0)) 



A * 2 ^ ^ (l+ 5 (0))(l+ 5 (0) + 2.g(V2/2))' 
which are plotted in Fig. 13. We see indeed that for 7 = 0.02 we have J-4 < J-3 and 
•7"4 < around Tq — 1/3. 

We summarize our results concerning regular patterns in Fig. ^|. We present the values 
of the Lyapunov function Tn as a function of 7 for N = 1, 2, 3, 4, 5, 6, 7, 8, and To = 1 /3, 
1, and 0. For T = 1 (Fig. |(c)) and ((Fig. g(d)), patterns of square symmetry (N = 2) 
have the lowest values of Tn for all values of 7 < 0.2, whereas the system favors patterns 
of different symmetries in different ranges of the value for 7 for Tq = 1/3. The second- 
harmonic resonance for r = 1/3, becomes less damped as 7 decreases, and thus the 
value of the self-interaction nonlinear coefficient g(0) becomes larger (see Fig. ||). In 
other words, the curve g(cji) has a wider flat center region and a sharper increase near 
cji = ±1, and therefore pattern structures with larger TV are favored. Table |l| shows the 
favored structures, and the corresponding ranges of 7. 

A couple of comments on the various patterns discussed above are in order, (i) A two- 
dimensional regular pattern structure with the spatial form C^j ex p(*% ■ x ) + c.c^j 
has N degrees of freedom. These N degrees of freedom appear as the phase of the complex 
amplitude Aj = A exp(9j) for j = 1, • • • , N. Among them, two correspond to spatial 
translations, whereas the other N — 2 degrees of freedom represents the phason modes 
(Golubitsky, Swift & Knobloch 1984; Malomed, Nepomnyashchii & Tribelskii (1989)). 
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For N = 3, the phase degeneracy for the single phason mode, which corresponds to 
hexagonal or triangular states, can be lifted by higher order nonlinear terms (Golubit- 
sky, Swift & Knobloch 1984; Muller 1993). Although it is beyond the scope of this article, 
it will be interesting to see how the phase relations of a spatial pattern is determined by 
higher order nonlinear terms, (ii) We predict that quasipatterns with large values of N 
occur at very small values of the linear damping coefficient 7. As discussed in the intro- 
duction, for very small values of 7 the mode quantization effect can be quite severe for 
finite-size systems. This implies that experimental verification of the quasipatterns with 
large values of N can be difficult f[| Obviously, a similar problem appears in numerical 
solutions of Faraday waves. 



The line ar stability of solutions (Eq. 3.41 ) of the standing wave amplitude equations 
(Eq. ( 3.38 )) can be determined by the spectrum of growth rates a of amplitude pertur- 
batio ns A j oc e at since all phase perturbations are neutrally stable as is obvious from 
Eq. ( 3.38 ). Equivalently, the linear stability of solutions can also be obtained by the 
eigenvalue spectrum of the matrix d 2 J-/{d\A 3 -\d\Ai\), linearized around t he s tationary 
solutions. In the context of the standing wave amplitude equations (Eq. 3.41), we can 
only consider very limited set of perturbations. The usefulness of the stability analysis 
is to identify unstable solutions. A stable solution in the context of the standing wave 
amplitude equations can, however, be unstable to other perturbations, such as a travel- 
ing wave mode or a perturbation with spatial variations. In the rest of this subsection, 
stable or unstable solutions are with respect to the perturbations that are allowed by the 
SWAE's. 

Let us consider a stationary solution with \ Aj | = ao given by Eq. (3.41) for j = 1, • • • , N 
and \Aj\ = for j = N + 1, • • • ,M (M > N), and small real perturbations bj to \Aj\ 
for j = 1, • • • ,M. On substituting the stationary solution with the perturbations into 
Eq. ( |3.38 ), we obtain the following linearized equations for by. 



1 



d T b 3 



1 



JV 

mi 1 3+ 9( c 



N 



bj - 2al Y 9{ c -ji)bi 



for j = 1, • • • , N, and 



7 



-d T bj 



1 




for j = N + 1,---,M. 

For parallel roll stationary solutions (N = 1), we have — 1 and, 

1 



7 



-drh = -26i, 



-dTbj = (1 - g(cij))bj 



for j = 2, • ■ ■ , M. Thus if g(cij) < 1 for some values of c\j, the parallel roll stationary 
solution is unstable. From the g(cji) curves in Figures [| ||, and [j| we conclude that 
parallel roll stationary solutions for Faraday waves in weakly dissipative fluids are not 
stable. 

The linear stability of square patterns (N = 2 and oq = 1/ y/l + g(0)) is determined 

f In laboratory experiments, mode quantization also depends on the nature of boundary 
conditions. Certain "soft" boundary conditions may relax the strict quantization requirements 
and allows access to the larger system regime than the actual size of the system (Douady 1990; 
Bechhoefer et al. 1995). 
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from the following linear system, 

l+g(0) 
2 7 

1 + 3(0) 
27 



9 T 6i = -61 - 5(0)62, 
d T b 2 = -6 2 ~ 3(0)6i, 



7 3 V 1 + 3(0) / J 



7 

for j = 3, • • • , M. Thus if -1 > g(Q) > 1 and 

, _ , gjcij) + g(c2j) 

A " _1 ~ 1 + 3(0) <0 ' 
square patterns would be stable. The first condition is satisfied as shown in Figures ||, 
[|, and||. Since |c2j| = Jl — Cy, we have 



3(ci J ) + 3(y / i-cf J ) 

A. 9 (cy) = l , 

where — 1 < C\j < 1. The growth rate A S(? of a standing wave perturbation in an arbitrary 
direction to the stationary square solution is always negative except for T = 1/3 and 
small values of damping coefficient 7 = 0.02. Interestingly, this is the parameter regime 
for To and 7 where hexagonal/triangular patterns are found to have lower values of the 
Lyapunov function than square patterns. The results of the stability analysis tell us 
that in this parameter range square patterns are in fact unstable, and thus correspond 
to a local maximum or a saddle point of the Lyapunov function. This result in turn 
guarantees that square patterns will not be seen in Faraday waves for this parameter 
range. The exact parameter range, however, cannot be determined from the stability 
analysis since we are only considering very restricted form of perturbations in the context 
of the standing wave amplitude equations. 

The linear stability of hexagonal or triangular standing wave patterns (N = 3 and 
do = 1/ yl + 23(1/2)) can be determined from the following linear system, 

1 + 2g(1/2 ld Tbl = - bl - 0(1/2X63 + 63), 



27 
23(1/2) 
27 

23(1/2) 
27 



a r 6 2 = -&2-3(l/2)(6 1 + &3), 
d T b a = -6 3 - 3(1/2) (61 + 6a), 



1 n h (1 9(cij) + g(c 2j ) + g{c 3j ) \ u 

for j — 4, • • • , M . Therefore, hexagonal or triangular patterns are unstable if 3(1/2) < 
-1/2, or g(l/2) > 1 or, 

, _-, a(cij) +3(c2j) + g(c 3 j) n 

Xht = 1 1 + 23072) ' 

for j — 4, ■ • • , M. Similar stability analyses can also be performed for quasipatterns 
(N > 4). 
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3.4. Envelope Equations 

We have assumed so far that Faraday wave patterns consist of a set of spatially uniform 
(i.e., no modulations) standing waves although we have considered the amplitude of the 
each standing wave to be slowly varying in time. It is conceivable, however, that the 
standing wave amplitudes may also have slow spatial variations, i.e, standing waves with 
spatially modulated amplitudes. In fact, the slow spatial variations of the amplitudes can 
be nicely incorporated into the amplitude equations within the framework first discussed 
by Newell fc Whitehead (1969)| for pattern formation in Rayleigh-Benard convection. 

Amplitude equations including slow spatial variations of the amplitudes are often called 
envelope equations. In principle one can also derive a set of coupled standing wave 
envelope equations for parametric surface waves from the hydrodynamic equations with 
the assumption that spatial variations of the amplitudes are of a slow scale of the order 
of e 1 / 2 , i.e., X = e 1 / 2 x and Aj — Aj(X.,T). With these assumptions, there will be no 
nonlinear terms involving spatial derivatives in the envelope equations up to C(e 3 / 2 ) fj] 
As a result, only linear terms involving spatial derivatives will appear in standing wave 
envelope equations (SWAE's), and the cubic nonl inear terms will be exactly the same 
as the standing wave amplitude equations (Eq. ( [3.34 )) (SWAE's). The linear terms 
involving spatial derivatives can be obtained by a similar perturbative expansion as for 
the SWAE's as well as from symmetry and invariance arguments, and the growth rate 
of the linearly unstable modes (Newell 1974; Cross & Hohcnbcrg 1993). We have chosen 
the latter way, which is much simpler. 

The possible scalar terms up to C(e 3 / 2 ) are (k 
Since the envelope equation for the standing wave amplitude Aj is expected to be invari- 
ant under the transformation kj — > —kj, the first term will not appear in the equation. 
The second term implies the same form of the transverse and longitudinal variations 
for a standing wave amplitude Aj. Since a standing wave pattern breaks the rotational 
symmetry of an isotropic two-dimensional system, the transverse and longitudinal spatial 
variations are qualitatively different, and thus the second term also cannot appear in the 

envelope equations |. Therefore, we only need to consider the last term (kj ■ Vxj Aj 
in the envelope equations, 



x 



Aj, V z x Aj, and 



A, 



7 dT 



A 3 +e,{k r ^x) 2 A 3 ^T) 



N 



9{c 3 i)\Ai 

!=i(i#i) 



A, 



(3.43) 



where £o is determined from the growth rate of the linearly unstable eigenmode as given 
byEq. (ph, 



a+(k) = - 7 + // 2 -p-l)(Go + 3r )/2] 2 . 
From the general arguments of Newell (1974) and |Cross fc Hohcnbcrg (1993) , we expand 



f Nonlinear terms involving spatial derivatives would appear in the envelope equations if the 
standing wave amplitude equations had quadratic nonlinear terms (Brandt 1989). 

X With the assumption of different slow spatial scales for the transverse and longitudinal 
modulations in analogy to Rayleigh-Benard convection, scalar terms of the following form are 

permissible in standing wave amplitude equations: 



kj-Vx - (i/2) (kj- ■ V x )' 



Aj with kj 



a unit vector perpendicular to kj. We do not consider such terms here since we are only 
computing the coherence length £o in the longitudinal direction. T he more general case of 
transverse modulations has been explicitly addressed by Zhang (1994) . 
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the growth rate u+(k) around the critical wavenumber k = 1 as 

1 , 



7-0 



Hence we hnd that r = I/7 and 

r / d 2 a + (k) 



dk 2 



tq ( Go 

V 



k=l 



3r 



The envelope equations now can be written as, 

2 



1^ 

7 dT 



A, 



Go + 3ro 
2V27 



A, 



N 



I A, 



Yl 9{cji)\Ai 



(3.44) 



Aj, (3.45) 



for j = 1, • • • , N. We have set / = 7 in the spatial derivative term since the difference is 
of higher order. 

The envelope equations (Eq. ( [3.45 )) are also of gradient form --gf- = —SC/5A*, with 
a Lyapunov functional 



r/x 



A' 



Go + 3Fq 



2\/2 7 



A, 



1 w 



JV 



£ 9{cji)\Ai 



(3.46) 



An important point about the standing wave envelope equations (Eq. 3.45) and the 
Lyapunov functional (Eq. |3.46| ) is that the coefficient of the spatial derivative term is 
much greater than one, i.e., 

since 7 <C 1 for weakly dissipative surface waves^ Physically £0 is a measure of the 

coherence length of the wave pattern. Thus Eq. ( 3.47 ) implies standing wave patterns 
of Faraday waves in weakly dissipative fluids have a coherence length much longer than 
their wavelength. This should be compared to Rayleigh-Benard convection, where we 
have £0 ~ 1 (Newell & Whitehead 1969; Cross & Hohenberg 1993), and to directional 
solidification of a binary alloy, where £ ^ 1 due to the extremely flat neutral stability 
curve (Mullins & Sekerka 1964). Because £ ^> 1, highly ordered square patterns as 
large as 30 — 40 wavelengths are often observed in Faraday waves in weakly viscous fluids 
(Gaponov-Gerkhov & Rabinovich 1990). 



4. Summary and Discussion 

We have studied in this chapter pattern formation in weakly damped Faraday waves by 
deriving a set of linear damping quasi-potential equations and by performing a multiscale 
asymptotic expansion close to onset. Standing wave equations have been derived that 
explicitly incorporate higher harmonic terms in the linear neutral solutions. These terms 
are seen to be important for the saturation of the wave amplitude (amplitude-limiting 
effect by the driving force). We have also studied in detail the effect of triad resonant 
interactions among capillary-gravity waves on nonlinear pattern selection. Triad resonant 
interactions are found to be the main reason for the appearance of square patterns in 
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capillarity dominated Faraday waves. By increasing the gravity wave component, the 
triad resonant condition is altered. As a result, square patterns can become unstable, 
and hexagonal or quasicrystalline patterns can be stabilized. 

The importance of higher harmonic terms in the linear neutral solutions also has 
implications for parametric surface waves in highly viscous fluids. Because of the larger 
threshold values of the driving force for highly viscous fluids, higher harmonic terms in 
the linear neutral solutions can be even more important than in the weakly damped case 
considered here. More than one higher harmonic term can be important in that case. 

Our results provide justification for the observed selected patterns of square symmetry 
near onset in fluids of low viscosity (Lang 1962; Ezerskii, Korotin & Rabinovich 1985; 
Tufhlaro, Ramshankar & Gollub 1989; Ciliberto, Douady & Fauve 1991; Bosch & van 
de Water 1993; Edwards & Fauve 1993; Miiller 1993). Without requiring any additional 
assumptions, Faraday waves close to onset are potential, and minimization of the associ- 
ated Lyapunov functional leads to square patterns in the capillarity dominated regime. 
This is in agreement with most experimen tal studies, with the exception of the work by 



Christiansen, Alstr0 m fc Levinsen (1995) who observed quasi-patterns where we predict 
squares. More recently, a systematic survey of pattern selection in a system of large as- 
pect ratio by Kudrolli & Gollub (1996) has again documented the transition to squares 
for capillary dominated waves, but also shown the transition to hexagons in the vicinity 
of r = 1/3, in agreement with our predictions. In the range in which we predict stable 
quasi-patterns, Kudrolli & Gollub (1996) observe hexagons instead. The origin of this 
discrepancy remains to be investigated. 

We are indebted to Maxi San Miguel for many useful discussions. This work was 
supported by the Microgravity Science and Application Division of the NASA under 
Contract No. NAG3-1284. This work was also supported in part by the Supercomputer 
Computations Research Institute, which is partially funded by the U.S. Department of 
Energy, contract No. DE-FC05-85ER25000. 
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Table 1. Patterns with the lowest values of the Lyapunov function and their corresponding 
ranges of 7. These values correspond to mixed capillary-gravity waves with To = 1/3. 



range of 7 favored pattern 



7 > 0.0809 square patterns 

0.0253 < 7 < 0.0809 hexagonal/triangular patterns 

0.0132 < 7 < 0.0253 eightfold quasipatterns 

0.0082 < 7 < 0.0132 tenfold quasipatterns 

0.0057 < 7 < 0.0082 twelvefold quasipatterns 

0.0042 < 7 < 0.0057 fourteenfold quasipatterns 




! A 




(b) 



Figure 1. Triad resonant interaction in parametric surface waves (a) 6^ as a function of IV 
(b) The mode B will resonate with the quadratic interaction of the standing waves Aj and Ai 
when 6ji =0$ ) . 




Figure 2. The self-interaction coefficient g(l) of the standing wave amplitude equations as a 
function of To with the linear damping coefficient 7 = 0.1 in (a), and 7 = 0.02 in (b). 
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Figure 3. The standing wave nonlinear coefficient g(cji) as a function of Cji for purely capillary 
waves (To = 1) with the linear damping coefficient 7 = 0.02 in (a), 7 = 0.05 in (b), 7 = O.f in 
(c), and 7 = 0.2 in (d). 
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Figure 4. The nonlinear coefficient g(cji) of the standing wave amplitude equation as a function 
of Cji for purely gravity waves (r = 0) with the linear damping coefficient 7 = 0.02 in (a), 
7 = 0.05 in (b), 7 = 0.1 in (c), and 7 = 0.2 in (d). 
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Figure 5. The nonlinear coefficient g(cji) of the standing wave amplitude equation as a function 
of Cji for gravity-capillary waves of To = 1/3 with the linear damping coefficient 7 = 0.02 in (a), 
7 = 0.05 in (b), 7 = 0.1 in (c), and 7 = 0.2 in (d). 
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Figure 7. The difference A32 of the values of the Lyapunov function for hexagonal/triangular 
patterns Tz and square patterns T2 as a function of To for 7 = 0.1 in (a) and 7 = 0.02 in (b). 




Figure 8. The differences in the values of the Lyapunov function for eightfold quasipatterns and 
hexagonal/triangular patterns, A43, in (a), and for eightfold quasipatterns and square patterns, 
A42, in (b). In both cases, 7 = 0.02. 




Figure 9. The values of the Lyapunov function Tn (N = 1, 2, 3, 4, 5, 6, 7, 8) as a function of 7 
for To = 1/3 in (a), To = 1 (capillary waves) in (c), and r = (gravity waves) in (d). The 
portion of (a) with small values of 7 is shown in (b). 



